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TWO-DIMENSIONAL SYMMETRICAL SUPERSONIC NO ZZLES 


By J. C. Eward and Lawrence R. Marcus 


SUMMARY 

Auxiliary Boundary conditions are derived to assure continuity of 
wall curvature in applying the method of characteristics to the design 
of two -dimensional symmetrical supersonic nozzles. An illustrative 
example is included. 


INTRODUCTION 

The method of characteristics was first applied by Prandtl and 
Busemann to the evaluation of supersonic flows in 1929 (reference l) . 
Since then, a number of expositions (for example, see references 2, 3, 
and 4) have described various applications of characteristics to the 
design of two-dimensional supersonic nozzles. The problem considered 
was to determine wall contours that would transform a uniform or source 
flow usually at Mach number of unity to a uniform shock-free flow at 
some higher Mach number. 

The wall contours so determined usually contain discontinuities 
in curvature which are tolerable for fixed geometry designs. In varia- 
ble Mach number flexible -walled nozzles, however, a discontinuity in 
curvature implies a multivalued stress at a point, which in turn can 
only occur if there is a restraining moment. The design of flexible- 
walled supersonic wind tunnels usually does not include provision for 
such a. restraining moment and hence the wall fairs the curvature across 
the discontinuity. If the characteristic diagram for the tunnel 
requires a discontinuity, then imperfections in the actual flow must 
appear . 

Even if a rigid-walled nozzle is used, the probability that the 
theoretical gas flow will occur is slight. Changes in the ratio of 
specific heats , or any condensation of the gas components, or imper- 
fections in the corrections for boundary layer can shift the character- 
istic pattern along the nozzle wall. The seriousness of such a s hif t 
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is minimized if the curvature is continuous and small. Here again noz- 
zles with inherent continuity of curvature from the characteristic 
pattern are desired. 

A plot of wall angles for a conventional supersonic nozzle including 
a discontinuity of curvature at the point of inflection is presented in 
figure 1. Also shown is a schematic plot on which the discontinuity in 
curvature has been eliminated. The auxiliary restrictions which must 
he Imposed in the design of two-dimensional supersonic nozzles to elim- 
inate discontinuity of wall curvature are presented herein. The work 
was completed during May, 1951, at the NACA Lewis laboratory. 

ANALYSIS AND DISCUSSION 

The problem is to determine the auxiliary conditions which must be 
satisfied to eliminate discontinuities of wall curvature in the design 
of supersonic nozzles. In conventional designs, a discontinuity occurs 
at the geometric point of inflection and also at the start of the test 
s ection . 

For convenience, two sets of characteristics i{r + and i{r“ will be 
utilized. In brief, 

* + + = v (1) 

is the Prandtl -Meyer expansion angle tabulated as a function of Mach 
number in reference 5. The quantity 

* + - r = o ■ ( 2 ) 

is the local flow deflection angle with respect to the axis . Either 
i[r + or ijr - characteristics represent Mach lines in the flow. The i{r + 
set (fig. 2) originates on the contour wall and moves from left to 
right as the characteristic is traced from the wall to the tunnel center 
line . The ifr - set is traced from left to right in projecting from the 
axis to the contour- streamline. Thus, at the wall, i{r + is always lar- 
ger than or . equal to if - . 

The characteristic point is that place on the wall where i{r + first 
assumes the test-section value. The point of inflection occurs at the 
position of maximum wa ll angle. 

The curvature along a streamline may be expressed mathematically 
■by 
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■where C is the curvature and s is the arc length. Continuity of 
curvature requires that 


d 0 _ _ df 

ds ds ds 


(2a) 


he continuous . 

At the characteristic point and downstream therefrom, i{r + is con- 
stant at the test- section value. . The curvature is then 

d0 _ dijr , 

ds ~ ‘ ^ ^ 

t 

downstream of the characteristic point. The point of inflection, on 
the other hand, has zero curvature so that at the point of inflection 

dfr + _ dfr~ 

ds ~ ds 


Thus, if the characteristic point coincided with the point of inflection, 
continuity of curvature would require that the local values of 

= - — = 0. A zero value of dijr - /ds at the point of inflection 

results only if di[r + /ds vanishes at some point on the wall farther upstream. 
Hence, at that upstream point, dfi/ds = -dijf-/ds. Because dif/ds must he 
positive to avoid compression waves, d©/ ds would then he negative and the 
nozzle would contain at least two inflection points. Rather odd-shaped 
walls therefore result which are undesirable for flexible nozzles. Thus 
the characteristic point should not coincide with the point of inflection 
if the wall curvature is to he continuous . (This coincidence does occur 
on most conventional nozzle designs.) Likewise, by similar arguments, the 
point of inflection should not he downstream of the characteristic point. 


The remaining possibility is for the point of inflection to he 
upstream of the characteristic point. No difficulty is experienced in 
achieving continuity of curvature at the point of inflection so long as 
di|r + diir“ 

Ijjj— and are continuous functions. Continuity can he achieved at 

the characteristic point also if ^ approaches zero from the upstream 


side. The curvature would then blend smoothly with its downstream value 
given by equation (3) . Thus t{t + is tangent at the characteristic point 
to its final value. 


Vanishing of at the characteristic point infers a vanishing 

— at the start of the test section (fig. 2) . Hence the discontinuity 

of curvature at the test-section wall Junction is likewise eliminated. 
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Thus the auxiliary conditions to assure continuity of wall curva- 
ture are: 

(l) The geometric point of inflection = Oj 

of the characteristic point. 


must he upstream 


(2) The derivatives 


d\jr 


and 


di|r _ 


must he continuous functions. 


ds as 

and. as a consequence i{r + must he tangent to its ultimate test-section 

value at the characteristic point. The continuity of ^ and 

ds dB 

also infers continuity of Mach number gradient along the axis . 


Nozzles with continuous wall curvature will always he longer than the 
shortest possible nozzle for a given Mach number and test-section size . 


An example follows for which the wall curvature is continuous . 

Plots of i(r + , 0, and Mach number distributions for an M = 4 nozzle 

are presented in figure 3. The flow was expanded from a uniform Mach 
number of unity at the throat in the usual manner (reference 2) to the 
wall point of inflection, where the curves of i[r + and ifr - have the 
same slope. The i{r + curve was then faired from tangency at the point 
of inflection to the constant-teSt-section value in a horizontal tangent 
at the characteristic point. The remaining values of Tjr“ and 0 were 
determined by the characteristic plots for the nozzle. As a consequence 
of the i{r + fairing at the characteristic point, 0 and i(r~ approach 
horizontal tangents at the test- section junction as does the wall Mach 
number . 


In practice, nozzles of predetermined length are often desired. 
Precise estimation of the nozzle length for a given test-section size 
is difficult when only the characteristic point and final Mach number 
are specified. An alternative design procedure is therefore useful. 

The flow- in the vicinity of the nozzle throat can be obtained by charac- 
teristics in the usual manner . A curve of either ijr + or Mach number is 
thus obtained along the axis which properly matches the initial super- 
sonic flow to the sonic flow at the throat. From the design Mach number 
and test-section size, the upstream vertex of the test diamond is then 
located. Continuity of dt]r+/ds and d\)r“/ds along the nozzle contour 
implies continuity of dijr + /ds and also dM/ dx along the axis . Either 
the i{r + value or the Mach number distribution along the axis should then 
be faired in a smooth curve to a horizontal tangent at the upstream vertex 
of the test diamond. Working from the axis the remaining wall contours 
are determined by the characteristic plots to give a nozzle of the 
required length and size. 


CONCLUSIONS 

Auxiliary bo undar y conditions were derived to assure continuity of 
wall curvature in the application of the method of characteristics to 
two-dimensional symmetrical supersonic nozzles. These conditions are: 
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1. The point of inflection must he upstream of the characteristic 
point. 

2 . The first derivatives of the characteristic values and 

ds 

must he continuous functions . Hence i^ + must he tangent to its 

test-section value at the characteristic point. This latter restriction 
infers that the wall slope approaches a minimum or horizontal tangent 
(that is, zero curvature) at the test-section junction. Also the axial 
Mach number gradient must he continuous and approach zero at the upstream 
vertex of the test diamond. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 16, 1951 
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Figure 3. - Variation of flow parameters for a Mach number 4 two-dimensional 

supersonic nozzle . 
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